Many-body interactions in quasi-freestanding graphene 
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The Landau-Fermi liquid picture for quasiparticles assumes that charge carriers are dressed by 
many-body interactions, forming one of the fundamental theories of solids. Whether this picture 
still holds for a semimetal like graphene at the neutrality point, i.e., when the chemical potential 
coincides with the Dirac point energy, is one of the long-standing puzzles in this field. Here we 
present such a study in quasi-freestanding graphene by using high-resolution angle-resolved pho- 
toemission spectroscopy. We see the electron-electron and electron-phonon interactions go through 
substantial changes when the semimetallic regime is approached, including renormalizations due 
to strong electron-electron interactions with similarities to marginal Fermi liquid behavior. These 
findings set a new benchmark in our understanding of many-body physics in graphene and a variety 
of novel materials with Dirac fermions. 



The many-body interactions in graphene at the charge 
neutrality point, where the Fermi surface of graphene is 
point like, are expected to differ from those of an ordinary 
metal. Due to the long-range nature of the electron- 
electron interaction, logarithmic velocity renormaliza- 
tions and a linear imaginary self-energy (ImS) are ex- 
pected by theoretical calculations, which differ from the 
case of a normal Fermi liquid [17 . The electron-phonon 
interaction also differs from that of a normal metal at 
the charge neutrality point, in that a finite self-energy 
is expected even in the absence of the ordinary electron- 
phonon "kink" [8l|9]. But despite the recent flurry of the- 
oretical works, there has been no experimental demon- 
stration of the self-energy in graphene at the neutrality 
point. 

Angle resolved photoemission spectroscopy (ARPES) 
is the ideal tool to probe electron dynamics and many- 
body interaction in graphene [10 as it can directly access 
the quasiparticle self-energy. So far however, most exper- 
imental studies of self-energy are on metallic graphene, 
i.e., when the chemical potential is far away from the 
neutrality point. In this regime many-body interactions 
are substantially screened [11 and can still be overall 
described within the Fermi liquid picture. Although 
these samples can be successfully hole doped to bring 
the Fermi level closer to the Dirac point, due to the ex- 
treme sensitivity of graphene electronic structure to dop- 
ing and/or disorder [12 , the resulting spectra are much 
broader than the as-grown samples [13 -15 making it dif- 
ficult to extract many-body information. Recently it has 
been shown that a freestanding monolayer graphene film 
can be found on top of a bulk graphite system, and that 
the monolayer is decoupled and suspended over the en- 
tire substrate [16]. Moreover it has been shown that 
its properties are nearly identical to a fully suspended 
graphene sheet [16]. An almost identical system can be 
found when graphene is grown on the carbon face of SiC, 
SiC(0001),p^. Due to the differing rotational orienta- 



tions of the adjacent graphene layers with respect to 
the substrate, layers can behave like isolated graphene 
sheets [nHE]. 

Here we present a systematic study of the self-energy 
of quasifreestanding graphene, grown on the carbon face 
of SiC. By comparing the experimental dispersions of 
quasifreestanding samples to theoretical calculations and 
highly-doped samples, we show that both the electron- 
electron and electron-phonon interactions go through 
substantial changes when the semimetallic regime is ap- 
proached. The electron-electron interaction exhibits log- 
arithmic renormalizations in the band dispersions, chang- 
ing the overall curvature of the band dispersions from the 
non-interacting or highly doped cases, and allowing us to 
extract the effective dielectric constant or fine structure 
constant a; while the imaginary self-energy is roughly 
linear in energy. We also show that the features of the 
electron-phonon interaction diminish significantly as the 
neutrality point is approached. 



RESULTS AND DISCUSSION 

Fig. lA shows slices of ARPES intensity through the 
Dirac cone of monolayer graphene, at the Brillouin zone 
corner, K. This conical dispersion with a point like Fermi 
surface is the typical feature of Dirac fermions in undoped 
monolayer graphene [191 HO] • The monolayer nature of 
this sample is also supported by the single vertical kz vs 
ky dispersion shown in panel b. If AB-stacked graphene 
were present, this dispersion would have a number of 
bands determined by the number of AB-stacked multi- 
layers present, with band intensity that oscillates as a 
function of kz [21 . To extract quantitative information 
about the doping level, in panel c we show the dispersion 
obtained by fitting the peak positions of the momentum 
distribution curves (MDCs), the intensity at constant en- 
ergy as a function of momentum. From the separation 
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between the tt bands at we conclude that the as- 
grown sample is slightly hole-doped with n = 8 x 10^^ 
cm~^, which is comparable in both sign and magnitude 
to that of exfoliated graphene and further supports 
the decoupled nature of this sample. From the slope of 
the fitted dispersion we can also extract the Fermi ve- 
locity, vf = 1.10 X lO^m/s, which matches the value 
measured by scanning tunneling microscopy [18^ . 

A close inspection of the electronic dispersion given 
in Fig. 2 A (red solid line), shows that there is a sig- 
nificant departure from a linear bare band (dotted red 
line), most noticeable in the vicinity of E-Ei? = 800 meV. 
This curvature can also be seen in the inset of panel a, 
where the ReS, defined as the difference between the 
bare band and the experimental peak position, is shown 
in the vicinity of the anomaly. This anomalous curva- 
ture is most likely due to a many-body effect, such as 
the electron-electron interaction, rather than a bare band 
effect. Indeed, both the LDA band curvature (dashed 
green line in the same panel) for undoped graphene and 
the valence band of highly doped graphene [TT] exhibit 
the opposite curvature to the experimental valence band 
reported here. Among the many-body interactions, the 
electron-electron interaction is the most likely candidate: 
the electron-phonon coupling should affect the dispersion 
at a much smaller energy (the phonon energy) and should 
have the opposite curvature as we observe here [H 111 [22] 
(see Fig. 4), with decreasing velocity at the Fermi level; 
similarly, the electron-plasmon interaction is negligible 
for neutral graphene, and should also have the opposite 
curvature in any case [5]. The electron-electron interac- 
tion therefore appears to be the source for the anomaly. 

This can be further confirmed by comparison between 
the experimental band and ab initio calculations for free- 
standing graphene, performed within the GW approxi- 
mation [23] (see supplementary information), also shown 
in Fig. 2A. The curvature of the electronic bandstructure 
is very similar when the electron-electron interaction is 
taken into account, as also seen in Ref. [23 . In con- 
trast, Ref. considers only the electron-phonon inter- 
action and does not predict this effect; as stated above, 
the electron-phonon interaction causes the dispersion to 
curve in the opposite direction, with decreasing veloc- 
ity at the Fermi level. In summary Fig. 2 reveals our 
primary result: that the origin of this finite curvature 
is due to the electron-electron interaction and results in 
an increase of the band velocity as the Dirac point is 
approached. 

Another interesting observation that follows from the 
comparison between the experimental and the theoreti- 
cal dispersion in panel 2a is the decrease of the velocity, 
i.e., the slope of the dispersion, with increasing dielec- 
tric screening. The LDA band, corresponding to infi- 
nite dielectric screening, has the smallest band velocity, 
while the ah initio data corresponds to a fully suspended 
graphene sheet and has the largest velocity. The experi- 



TABLE I: Summary of band velocities (m/s) 

Undoped Highly Doped 

Dispersion YFermi ^High ^Fermi ^ High 

LDA 0.85x10^ 0.89x10^ 0.51x10^ 0.61x10^ 

Experiment 1.10x10^ 1.13x10^ 0.49x10^ 0.77x10^ 

Ah Initio GW plus el-ph 1.48x10^ 1.80x10^ 0.58x10^ 0.74x10^ 

^High is the band velocity at approximately 300meV, 
slightly higher binding energy than the phonon energy. 

mental band velocity lies between these two extremes due 
to the finite dielectric screening of the graphitic substrate 
(the finite doping of the experimental data may also play 
a role). Values of the band velocities near the Fermi level 
are listed on the left side of table 1. 

In contrast to undoped graphene (panel a) the effect of 
the electron-electron interaction and dielectric screening 
is reduced for highly doped graphene (panel b), resulting 
in a better overall agreement with the ab intio calcu- 
lations (a summary of band velocities can be found in 
table 1). The extent to which this might be true can 
be examined by considering the inverse screening length 
for doped graphene, which is given by qs = Aakp [24J. 
Our as-grown sample has a Fermi vector of /c^ = 0.005 
A~^, while the highly-doped sample has a Fermi vector 
of kp = 0.23 A~^. Using our experimentally determined 
value oi a = 0.4 (discussed below), we find length scales 
corresponding to 790 A and 17 A for the as-grown and 
highly-doped samples, respectively. In addition to veri- 
fying the validity of the ah initio calculations, a close in- 
spection of the highly-doped dispersions shows that they 
curve in the opposite direction as the as-grown sample. 
These results clearly indicate that the many-body physics 
of freestanding graphene near the charge-neutrality point 
is dramatically different from highly-doped graphene. 

To further investigate the implications of these find- 
ings, in Fig. 3a we consider Hartree-Fock and renor- 
malization group corrections to the self-energy to ac- 
count for the high energy band curvature observed in 
Fig. 2a. It was proposed that, in contrast to the lin- 
ear real self-energy of a Fermi liquid where quasiparticles 
decay via intraband transitions, when graphene is near 
the charge neutrality point, the electron-electron interac- 
tion leads to a logarithmic correction of the self-energy 
[TJ [25] . This logarithmic self-energy can be approximated 
as Rei;ee(^) = I^fI^ ~ ln(|^^3^|), where kp is the 
Fermi wavenumber along the T — K direction, ko = 1.703 
A~^ is the Dirac point wavenumber, \kD — /^f| = 0.005 
A-\ kc = 1.703 A-i is a cutoff, and the fit param- 
eter a is the dimensionless fine structure constant of 
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graphene, a = Consequently, the dispersion rela- 

tion is modified as E{k) = —Vp\k — /c^l + ^ee{k). We 
obtain a bare band velocity of Vp = (0.86 ± 0.02) x 10^ 
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FIG. 1: The band structure of monolayer graphene. (a) ARPES spectra of monolayer graphene, showing several slices through 
the Dirac cone of monolayer graphene. (b) The /c^ -dependence of photoemission intensity shows a single line, signature of 
monolayer graphene. (c) The red solid line shows the peak positions of the MDC fits, and the dotted red line represents the 
linear fit of the dispersion. The position of the Dirac point is slightly above the Fermi level. 




FIG. 2: Self-energy effect of the electron-electron interaction, (a) Extracted dispersions from as-grown C-face graphene (red), 
and from ab initio GW plus el-ph calculations including el-ph and el-el interactions in suspended graphene (black). Straight 
red and black dashed lines are given as guides to the eye, which intersect the dispersions at -1.2 eV and -0.2 eV, to show how 
the experimental dispersions (solid lines) deviate from linearity (straight dotted lines). The LDA (local density approximation) 
band is given as a green dashed line. Inset: Deviation of the experimental dispersion from linearity (red dotted line subtracted 
from the solid red experimental dispersion), (b) Extracted dispersions for highly doped graphene (blue). Doping is obtained 
by K-deposition and the resulting charge density is more than three orders of magnitude higher than undoped graphene. Black 
solid line shows the ab initio GW plus el-ph calculated dispersion. The LDA band is shown by the green dashed line. 



m/s (which is smaller than the experimental Fermi veloc- 
ity of 1.10x10^ m/s) and an effective dielectric constant 
e = 6.4 ± 0.1 (or a = 0.40 ± 0.01). The two param- 
eter fit was performed over the energy range [-1.3 eV,- 
0.4 eV] to avoid contributions from the electron-phonon 
kink, in addition to the constraint that the band passes 
through E{k = kp) = 0, to account for the low bind- 
ing energy deviations in the vicinity of the phonon kink 
(which will be discussed later), and some finite-resolution 
effects near the Fermi level. The very good agreement 
between the logarithmic fit and the measured dispersion 
(panel 3a) points to a clear signature of renormalizations 
due to electronic correlations. This self-energy takes the 
same form as that of a marginal Fermi liquid, although 
higher order corrections might lead to small deviation 
from this picture [24 . Note that, although obscured by 
the presence of the phonon kink, electron-electron renor- 
malizations remain apparent at low binding energy: the 
experimental band velocity remains larger than that of 
the LDA band, and the electron-electron interaction is 
the only interaction known to result in an increase in the 
Fermi velocity. 



The extracted value for the dielectric screening falls 
within the large range of values reported in the liter- 
ature, which are influenced by the choice of substrate. 
These values range from e = 15 for a freestanding 
layer [26], to e = 4.4 for monolayer graphene on SiC, 
although this value increases greatly in the presence 
of an inter facial carbon layer [2 7 . Theoretical calcula- 
tions typically discuss a dielectric constant of the form 

^ — (^^substrate H~ ^vacuum) / — (^^substrate H~ l)/^? with 

an additional term for screening from the graphene film 
itself: calculations expect e = 4 for fully suspended 
graphene [28^. A comparison between our data and these 
calculations lead us to expect that the graphitic substrate 
has a dielectric constant of roughly e = 7. This value is 
of the same order of magnitude as previous studies, al- 
though the real dielectric constant of graphite is a source 
of controversy with measured values ranging from 2 to 
13 [29 . 

The presence of a logarithmic term in the self-energy, 
and hence the departure from the Fermi liquid picture, 
should also be manifested in the ImS, predicted to scale 
linearly with energy [2 . This can be observed in Fig. 
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FIG. 3: Comparison of the experimental dispersions to Hartree-Fock and renormalization group self-energy, (a) The experi- 
mental dispersion (red) and theoretical fit (black) given by the electron-electron self-energy correction to the bare band. The 
fit is given by E{ky) = 5.63(A:^ - kp) - 0.563{ky - /cf) ln(| |). (b) ImE extracted from Voigt fits to MDCs. The dotted 

line is a guide to the eye. 



3b, where a linear ImXl is observed over a large energy 
window. The ImS is extracted from the experimental 
data using the standard fitting procedure by Voigt line- 
shapes of the raw MDCs [30], where the Voigt func- 
tions are defined as the convolution between a Gaus- 
sian, to account for the resolution broadening, and a 
Lorentzian, the expected MDC lineshapes in the presence 
of a /c-independent spectral function. This fitting proce- 
dure allows one to isolate the intrinsic contribution to 
the self-energy, by separating extrinsic parameters such 
as experimental resolution and sample inhomogeneities 
[30 . Note that the linear self-energy extends up to more 
than 1 eV, as expected in the case of conical dispersion 
[31], where the energy range of validity of the logarith- 
mic renormalizations roughly coincides with the energy 
where the band structure loses its linear dispersion rela- 
tion [32]. However, we should also note that the electron- 
phonon interaction adds a significant contribution to the 
self-energy, which in undoped graphene is predicted to 
increase in binding energy in a similar manner [9 . 

These findings are in line with infrared spectroscopy 
studies, where electronic correlation effects were reported 
to modify the Fermi velocity [33 . It also supports the 
view that electronic compressibility measurements [34] 
can be explained in terms of electronic interactions, as 
discussed by Sheehy and Schmalian [35] . 

In Fig. 4 we focus on the low energy regime, in the 
range of 150-200 meV from E^?, where another deviation 
of the dispersion relation from linearity is observed. This 
deviation, also known as kink, has been widely studied 
for highly doped graphene and is a signature of electron- 
phonon interaction [S] (22. What is surprising about 
these data is that, while for doped graphene the effect 
of the electron-phonon interaction is strongly manifested 
in the dispersion (panel b) and in both the ReS and ImS 
(panels c and d respectively) in the form of a peak and a 
step at the phonon energy respectively, it is strongly re- 
duced for neutral graphene, where only a tiny deviation 
from linearity can be discerned (panel a and c). This 
translates into a featureless ImS (panel d). 

This strong reduction of the kink for neutral graphene 
is well reproduced by theoretical calculations of the 



electron-phonon self- energy [9] where electron-electron in- 
teractions (and others) are not included (panels c and 
d). On the other hand, the slopes of the theoretical 
and experimental self-energies differ greatly when only 
the electron-phonon self-energy is considered, since the 
electron-phonon interaction is only one component of the 
total self-energy. It is important to note that the absence 
of a feature in the self-energy does not necessarily imply 
that the self-energy is zero. Since each of these inter- 
actions have self-energy components that increase with 
binding energy, we can only compare the distinct features 
of these interactions. In the case of the electron-phonon 
interaction, undoped graphene has a finite self-energy, 
but no step in the ImS at the phonon energy, while the 
kink in the ReS is diminished as well lOj. These effects 
are clearly indicated in our data, and differ from the be- 
havior of ordinary metals. 

CONCLUSION 

In summary, we have demonstrated that many-body 
interactions in monolayer graphene are substantially 
modified when the neutrality point is approached, with a 
departure of the self-energy from ordinary metallic sys- 
tems. The electron-electron interaction results in novel 
renormalizations, while the effects of the electron-phonon 
interaction diminish as a function of decreasing doping. 
These findings are the result of the point-like Fermi sur- 
face and linear bands of monolayer graphene, which set 
the many-body interactions in graphene apart from those 
of ordinary metals. 

MATERIALS 

The graphene layers were grown in a closed rf induc- 
tion furnace at a temperature of 1550° C (see Ref. [36] 
for details). The sample was annealed to 1000° C in 
ultra-high vacuum (UHV) prior to measurement. Our 
ARPES investigation was performed at beamline 12.0.1.1 
at the Advanced Light Source, at a pressure better than 
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FIG. 4: Doping dependence of the electron-phonon interaction. The dashed hnes are guides to the eye. (a) Dispersion (sohd 
red Une) extracted along the F-K direction for the as-grown graphene near zero doping. A small deviation from linearity (light 
blue dashed line) is indicated by the arrow, (b) Dispersion (solid blue line) extracted from the F-K direction of highly-doped 
graphene. A kink due to electron-phonon coupling is indicated by the horizontal arrow (the light blue dashed line is a guide 
to the eye), (c) Comparison of ReS for doped graphene (blue), undoped (red) graphene, and el-ph calculation (orange). To 
approximate the bare band, a linear offset has been subtracted, forcing each curve to intersect zero energy at -0.4 eV and 
eV. (d) Comparison of ImS for doped graphene (blue), undoped graphene (red), and electron-phonon calculation (orange). For 
ease of viewing, an offset was removed from the highly-doped data. 



3xl0~^^ torr, with a sample temperature of 15°K, pho- 
ton energies of 42 eV and 50 eV. All of the ARPES 
measurements in this study were performed on a sin- 
gle rotational domain of a single graphene sample. The 
lateral extent of this rotational domain is less than 30 
/imx80/im, the approximate dimensions of the ARPES 
beamspot. This sample was then electron-doped in situ 
by potassium deposition with an SAES Getters alkali 
metal dispenser, allowing us to study the same position 
on the sample for both the as-grown and highly-doped 
cases. Ah initio calculations including the contributions 
of electron-electron and electron-phonon interactions to 
the self-energy were performed as in the supplementary 
information. 
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